火星轨道交会自主导航与制导方法

对火星采样返回任务中的火星轨道交会自主导航和制导技术进行了研究。采用光学自主导航敏感器测量的火星中心方向和视半径，相对敏感器测量的相对位置等观测量，设计了导航滤波器同时估计轨返组合体和上升器的轨道。在导航滤波器设计中，针对光学自主导航敏感器更新频率远低于滤波解算频率的问题，设计了一种连续观测量构造算法，确保每个滤波周期均可进行测量更新，以提高导航精度。基于导航滤波器估计结果，采用T-H制导设计了4脉冲共椭圆交会策略实施轨道控制，从而构成近程交会自主导航和制导方案用于完成火星轨道交会任务。通过数学仿真校验了所提出方法的有效性。     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

同步双小行星系统共振轨道设计

研究了同步双小行星系统中共振轨道的设计方法及演化规律。首先，基于双椭球模型建立探测器运动方程，并给出共振轨道初值选取方法。然后，利用改进并行打靶法，提出一种双小行星系统平面共振轨道两步修正方法。同时结合稳定性理论及分岔理论，给出双小行星系统三维共振轨道生成和延拓方法；最后，以双小行星系统1999KW4为例，设计了共振比为1∶1,1∶2,1∶3,1∶4,2∶3的平面和空间共振轨道族，并分析了共振轨道的特性及轨道周期和轨道能量的变化规律。给出的双小行星系统中共振轨道的设计方法具有普适性，对未来双小行星系统探测任务中的轨道设计具有一定的参考意义与借鉴价值。     

基于TLE的弹道系数计算方法及应用分析

针对仅使用两行要素（Two Line Element，TLE）作为数据源的应用需求，研究了基于TLE轨道衰减的弹道系数计算方法。介绍了一种常用的基于两组TLE的直接计算法，分析TLE选取间隔对结果精度的影响；提出了一种基于多组TLE的迭代计算方法，以降低异常TLE对计算结果的影响；从弹道系数计算效果、在再入预报中的应用等方面对这两种方法进行比较分析。结果表明，两种方法各有优劣，基于多组TLE的迭代计算法稳定性更高、受TLE精度的影响更小；由于数据区间更短，基于两组TLE的计算结果对短期轨道衰减特性反应得更准确，用于临近再入时的预报效果更好。     

地球静止轨道微波辐射计技术

静止轨道微波辐射计可实现全天时、全天候、高频次云雨大气观测,在台风、流域性强降水探测和预报方面发挥重要作用。但在波束扫描、多频复用、系统定标等方面的运用存在较大困难,世界上尚无在轨应用先例。现从探测需求出发,对静止轨道微波辐射计进行了系统方案设计。提出了卫星平台与辐射计部件快慢结合扫描的方法,以解决波束扫描和系统定标难题;提出了低频选入射角层叠式布局准光学频段分离方法,以解决多频段复用难题。研制了微波辐射计原理样机,通过实验室和外场试验验证了系统设计,为静止轨道微波探测卫星研制提供了一定指导。     

超低轨道卫星应用离子电推进技术方案

低地球轨道大气环境对诸如科学探测和对地观测卫星的阻尼作用十分明显，而且阻尼随太阳和地磁活动以及昼夜、季节交替变化范围宽。为了保证卫星轨道精度或飞行状态满足任务要求，需要利用推进系统对卫星受到的阻尼进行实时或间歇式补偿以实现轨道或飞行状态的保持。针对轨道高度220~268km的无拖曳飞行和轨道维持应用，基于卫星轨道阻尼变化和有效载荷指标要求分析，研究确定了离子电推进技术指标、推力调节方案、系统组成、推力控制方案和在轨应用策略，并对推力调节方案进行了试验验证。结果表明，与无拖曳飞行卫星任务匹配的离子电推进指标为推力调节范围1~20mN，推力分辨率优于12μN，与对地观测卫星轨道维持任务匹配的指标为推力调节范围1～25mN，推力分辨率100μN。研究提出的针对超低轨道卫星应用需求的高精度推力连续调节离子电推进技术方案，具有工程任务针对性和参考价值。     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

Effects of in-track maneuver on Sun illumination conditions of near-circular low Earth orbits

By developing approximate analytical models considering the J₂ perturbation, the effects of an in-track maneuver on the orbital Sun illumination conditions of near-circular low Earth orbits are analyzed. First, two approximate models for the variations in orbital sunshine angles are developed, one for variations at a given time and the other for variations at a given argument of latitude. Next, two approximate models for variations in orbital arc in Earth shadow are developed, one considers the small eccentricity and the other uses the zero eccentricity. Finally, the developed approximate models are applied to analyzing the Sun illumination conditions of a typical in-track maneuver mission on a near-circular low Earth orbit. From the results obtained, three major conclusions can be drawn. First, the variations in orbital sunshine angles at a given time may reach tens of degrees when the drifting time reaches hundreds of orbital periods, and the approximate model for that situation cannot effectively approach the numerical results. Second, the variations in orbital sunshine angles for any given argument of latitude are only a couple of degrees even when the drifting time reaches 500 orbital periods, and the approximation model developed can effectively approach the numerical results. Third, for variations in orbital arc in Earth shadow, the approximate model considering the small eccentricity has simple expressions and can effectively approach the numerical results; in contrast, the approximate model using the zero eccentricity has relatively worse precision.     

30cm离子推力器在轨环境下的工作点及热控措施研究

栅极间距变化是影响离子推力器在轨环境下从冷态条件正常点火启动的重要因素,同时也决定了离子推力器的在轨工作时机和热控实施策略。本文采用有限元仿真与地面热平衡试验验证相结合的方法,建立起30cm离子推力器有限元分析模型并进行了模型校验,之后对离子推力器在轨受太阳光照影响的栅极温度场分布和间距变化,以及推力器在5kW工况下的三个典型温度点所对应的栅极间距变化进行了仿真分析,最后考虑了主动热控干预对推力器最恶劣工作点的栅极间距变化影响。结果显示：纯太阳光照影响下的栅极组件存在周期性温度变化,栅极最大温差可达到100℃,栅间距缩小量在0.06mm~0.16mm范围内波动;在太阳光照基础上实施60W的主动热控后,栅极最大温差降低至60℃,栅间距缩小量波动范围则变为0~0.03mm;栅极最高温度点和最低温度点分别是推力器冷态启动最容易和最困难的两个工作时机点,两点所对应的启动后屏栅和加速栅最小间距分别为0.22mm和0.04mm;在10W、70W和120W的热控加热功率下,从最低温度点启动后的屏栅和加速栅最小间距分别为0.06mm、0.20mm和0.29mm;采取主动热控措施能够有效降低推力器工作过程中的栅极热形变位移峰值,且加热功率为120W即温控点温度为50℃的主动热控可以满足30cm离子推力器在轨冷态启动时的0.25mm安全栅极间距要求。     

深空测控网干涉测量系统在“鹊桥”任务中的应用分析

在"嫦娥4号"任务的第一阶段—"鹊桥"阶段,北京航天飞行控制中心利用佳木斯及喀什深空站对"鹊桥"进行了干涉测量观测,获取了实时与事后的高精度测角观测量,有效支持了任务的实施。两深空站需同时完成测控任务,无法交替射电源观测来进行系统差标校,基于此系统采用了长时间隔、在航天器观测前及双站结束后观测射电源的标校方法,在地月转移段、月球至L2转移段、Halo轨道形成段开展了多次干涉测量观测,所获得的时延、时延率结果直接应用于事后联合轨道确定,结果表明:深空网的时延观测精度约为3 ns。